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Abstract 

The formation and spectrum of two-particle bound states are investigated within a simple 
relativistic quantum field model with the Yukawa-type interaction. Within this approach, 
relatively weakly interacting quarks and gluons form stable bound states under the analytic 
confinement. By introducing a minimal set of free parameters (the quark masses, the coupling 
constant and the confinement ’radius’) and by solving the Bethe-Salpeter equation with one- 
gluon exchange, we satisfactorily explain the experimental data for most known mesons and 
the glueball. 


1 Introduction 

The conventional theoretical description of colorless hadrons within the QCD implies 
that they are bound states of quarks and gluons considered under the color confinement. 
Another realisation of the confinement is developed in [1] by using an assumption that 
QCD vacuum is realized by the self-dual homogeneous vacuum gluon field that can 
lead to the quark and gluon confinement. Hereby, propagators of quarks and gluons 
are described by entire analytic functions in the p 2 -complex plane [2], i.e. the analytic 
confinement takes place, although there exists a prejudice to the analytic confinement 
conception (e.g., [3]). 

In our earlier paper [4] we considered simple scalar-field models to clarify the 
pure role of the AC in formation of the hadron bound states. In doing so, we have 
demonstrated just a mathematical sketch of calculations of the mass spectrum for the 
’’mesons” consisting of two ’’scalar quarks” by omitting important quantum degrees of 
freedom such as the spin, color and flavor. Nevertheless, these models gave a quite rea¬ 
sonable sight to the underlying physical principles of the hadron formation mechanism 
and have resulted in qualitative descriptions of the spectra of ” mesons”, their Regge 
trajectories and ’’glueballs” [5]. 

It seems reasonable to consider below a more realistic extention of our approach by 
taking into account the color, flavor and spin degrees of freedom for constituent quarks 
and gluons. In doing so, we estimate the spectra of pseudoscalar and vector mesons as 
well as the glueball lowest state. 


2 Model 

Conventionally, considering this problem within QCD one deals with complicated and 
elaborate calculations, because the confinement is achieved as a result of strong inter- 
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action - by involving high-order corrections [6]. Then, there arises a problem of correct 
and effective summation of these contributions. 

On the other hand, the use of a QFT method is effective when the coupling is 
not large. Then, lower orders of a perturbative technique can result in satisfactory 
accurate estimates (e.g., in QCD) of observables. Particularly, one can effectively use 
the one-gluon exchange approximation. 

We consider a relativistic physics and for the hadronization processes use the Bethe- 
Salpeter equation, because the binding energy is not neglible and the relativistic cor¬ 
rections are considerable. In doing so, we use a minimal set of free parameters. 

The Lagrangian of the Yukawa-type interaction reads: 


c = 


Ja/3^/3 


+ 
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(i) 


where \I/“^(x) = Tg°j° r flavour ( x ) an d 0^0*0 = dvector( x ) are the quark and gluon fields. 

We guess that the matrix elements of hadron processes at large distance (in the 
confinement region) are integrated characteristics of the quark and gluon propagators 
and the interaction vertices. Therefore, taking into account the correct symmetry 
features is more important than the tiny details of these functions. 

Our aim is to suggest the most plain forms of these propagators which keep the 
essential properties and result in a qualitative and semi-quantative description of the 
hadron spectra. We consider the following propagators for quarks and gluons: 


0 J p z +m z 
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where A - the confinement energy scale, m - the quark mass for given flavour and 0 < 
c o(m) < 1 is a smoothly decreasing function. These are entire analytic functions and 
could serve reasonable approximations to the explicit propagators in the hadronization 
region. 

Consider the partition function and take an explicit integration over 0-variable: 
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JJ 5T5T e ~ CFiq,M J 5(j) 


II 
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( 3 ) 


where the two-quark and two-gluon bound states in the one-gluon exchange mode are 
generated by terms 


jC 2 
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3 Light and Heavy Mesons 

Below we shortly describe important steps of our approach for two-quark bound states. 
First, we allocate the one-gluon exchange between quark currents and go to the relative 
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co-ordinates in the centre-of-masses system. Then, make a Fierz transformation to get 
the colorless bilocal quark currents and introduce an orthonormal basis {Uq(x)} indiced 
by the quantum numbers Q = {nKlm}, where l stands for the orbital quantum number. 
We diagonalize C 2 on colorless quark currents and use a Gaussian representation by 
introducing auxialiary meson holds Bj y. Then, we apply the Hadronization Ansatz to 
identify Bj^(x) as meson helds with quantum numbers J\f = {JQff}. The partition 
function for mesons reads: 



DB+f e -\( B MW+9 2 T±V N SVMS)B M )+W I lB M } 


( 5 ) 


where S - the quark propagator and V)y = Y j j dy Uq{%))\J D(y) exp{f 0} is a vertice 
function. The interaction between mesons is described by Wj[B^-\ = 0[Bj/\. 

The diagonalization of the quadratic form is equivalent to solution of the Bethe- 
Salpeter equation (in the ladder approximation) on the orthonormal system f/y: 

g 2 Tr{Vj^SVM'S) = jj dxdy U JQ (x ) g 2 U p (x,y) Uj> Q >{y) = - 5 QQ ' 5 JJ ' X JQ (-p 2 ). (6) 

Then, the meson mass is dehned by equation: 


1 — X y(JWy) — 0. 


( 7 ) 


The polarization kernel fl p (x,y) is real and symmetric that allows one to use a 
variational technique. For simplicity we consider orbital excitations: n — 0, l > 0. By 
solving the variational equation for the meson mass 
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, PJ = 1 1 )of> Xj = {+1,-1} for J={P,V} 


we fit the known experimental data and the optimal set of parameters is obtained: 


a s = — = 0.324 , A = 668MeV, m u = rn ( i = 179MeV, 
47T 

m s = 202MeV, m c = 781MeV, mb = 4324MeV, 


( 8 ) 


3.1 Meson Ground States 

With these parameters first we estimate the pseudoscalar and vector meson masses for 
the ground state l — 0 in the wide range from light 7 t( 138) to heavy Y(9460) represented 
in Table 1. Note, our parameter a s is close to a particular value a s = 0.34 known from 
the r-lepton decay. In the present paper we consider fixed coupling constant in the 
range from hundred MeV to ~ lOGeV, although there exists another point of view to 
introduce a running a s dependent on the energy scale [8, 9]. 
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tt(140) 
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£>(1870) 

1929 

p(770) 

770 

K* (892) 

891 

p (547) 

547 

£,(1970) 

2012 

u( 782) 

782 

£*(2010) 

2001 

?7c(2979) 

3053 

£(5279) 

5356 

$(1019) 

1004 

£*(2112) 

2081 

Vb( 9300) 

9454 

£,(5370) 

5398 

J/T(3097) 

3097 

£*(5325) 

5356 

K (495) 

495 

£ c (6400 ± 400) 

6154 

T(9460) 

9460 




Table 1: Estimated masses (in MeV) for the pseudoscalar and vector mesons. 


3.2 Meson Excites States 


Compared to the ground states, the orbital excitations l > 0 take place in larger dis¬ 
tances and should be less sensitive to the short-range details of the quark and gluon 
propagators. A correct description of the mesonic Regge trajectories can serve an addi¬ 
tional testing ground for our model, for the approximated propagators. Our estimates 
(solid lines) for the p-meson and K -meson families excitations plotted versus J = 1 + S 
is given in Figure 1 in comparison with the experimental data (dots). 

As is expected, our model describes well the Regge trajectories of mesons. The 
slope of the Regge trajectories is very sensitive on the confinement scale parameter A 
and its optimal value is fixed near A = 668 MeV. 



Figure 1: Mj plotted versus the quantum number J = l + S (the Regge trajectories) 
for the p- and iW-meson families. Solid lines correspond to our estimates and dots 
depict the experimental data from Particle Data Group 2004. 


4 Glueball 


There may exist the glueballs, bound states of gluons without any constituent quarks. 
The experimental status of glueballs is not clear [10], although it is evident that the 
structure of QCD vacuum plays the main role in their formation. We suppose that 
the lowest-mass glueball may be approximated by the pure ’’gluon-gluon” bound state 
with scalar quantum numbers J pc = 0 ++ . 

Omitting details, the equation for the lowest-state glueball mass reads 


1 - fl {-Ml) = 0 
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With obtained values {aq = 0.324, A = 668 MeV} we estimate 


M g = 1433 MeV. (10) 

This result may be compared with other predictions of the glueball mass: 

M Gs = 1580 ± 200 MeV [Y.A.Simonov , 2000], 

M Gf = 1530 ± 200 MeV [H.Forkel , 2001]. 


Concluding we remark that our simple relativistic field model with an analytic 
confinement of quarks and gluons gives a quite reasonable sight to the underlying 
physical principles of the hadron formation mechanism. It allows one to describe 
qualitatively the meson spectra in the wide range from light 7 t( 138) to heavy Y(9460). 
By using only a minimal set of fundamental parameters (the constituent quark masses 
rtif , the coupling constant a s and the confinement scale A) we are able to describe 
the meson ground states, the Regge trajectories of orbital excitations and the low- 
mass glueball state with sufficient accuracy. Hereby, the coupling constant is relatively 
small that makes the use of one-gluon exchange modes reliable for the consideration. 

The author thanks Ts.Baatar, G.V.Efimov, S.B.Gerasimov, E.Klempt, S.N.Nedelko 
and W.Oelcrt for useful discussions and comments. 
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